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ABSTRACT 


In  the  first  phase  of  this  contract,  we  have  developed  an  integrated  approach  to  using  wav  on  minor 
parameters  and  their  covariance  properties  from  a  collection  of  local  arrays  for  estimating  location, 
along  with  an  uncertainty  ellipse.  Hypothetical  wavennnber  estimators  and  their  uncertainties  are 
used  as  input  to  a  Bayesian  nonlinear  regression  that  produces  fusion  ellipses  for  e\ent  locations 
using  probable  configurations  of  detecting  stations  in  the  proposed  global  infrasound  array.  The 
netw  ork  capability'  is  characterized  as  a  function  of  separate  local- array  characteristics,  including 
signal-to-noise  ratios,  bandwidth,  array  geometry,  local  correlation  and  coherent  interfering  signals. 
A  summary  map  displays  the  average  areas  of  the  90%  posterior  probability'  ellipses  for  each 
hypothetical  location,  assuming  a  random  configuration  of  detecting  stations. 

In  the  second  phase  of  the  project,  we  are  developing  local-array  parameters  that  will  be  used  as 
input  for  estimating  the  capabilities  of  the  global  International  Monitoring  System  (IMS).  A 
small  array  theory  has  been  given  in  previous  work  that  characterizes  the  detection  probabilities 
and  large  sample  variances  of  the  local- array  optimal  maximum  likelihood  detectors.  We  are 
w  or  king  on  assessing  the  local-arry  performance  of  the  multiple  signal  F-statistic  as  well  as  those 
of  alternative  high  resolution  detectors  produced  by  Capon  (1969)  and  the  multiple  signal 
classification  (MUSIC)  algorithm  proposed  by  Schmidt  (1979).  F  or  pure  and  mixed  infrasound 
signals  from  tw  o  explosions,  w  find  that  all  statistics  have  comparable  resolution.  The  F-statistic 
retains  a  number  of  theoretical  advan  tages,  namely  (1)  a  known  large  sample  distribution  that 
yields  detection  and  false  alarm  probabilities  (2)  direction  of  arrival  (DOA)  estimators  with  means 
and  covariance  matrix  determined  by  the  Cramer-Rao  lower  bound  and  (3)  easily'  estimated 
signal-to-noise  ratios.  Hence,  the  maximum  likelihood  procedure  produces  the  input  necessary  for 
evaluating  the  global  location  performance  of  the  IMS. 

With  the  above  in  mind,  current  efforts  are  focused  on  analyzing  events  in  the  infrasound  database 
located  at  the  Center  for  Monitoring  Research.  We  are  planning  to  develop  from  the  CMR  R&D 
T  est  Bed  Infrasound  Wveform  Library  information  on  detection  probabilities  and  estimated 
variances  for  real  IMS  type  arrays  as  a  function  of  bandwidth,  signal-to-noise  and  geometry.  These 
characteristics  are  needed  for  input  into  the  multiple  station  global  location  simulation. 

KEY  WORDS:  maximum  likelihood,  Bayes,  multiple  signal  estimation  and  detection,  nonlinear 
regression,  F-Statistics,  high  resolution,  MUSIC 
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OBJECTIVE 


Monitoring  explosive  events  using  a  collection  of  small  infrasound  arrays  can  lead  to  improved 
detection  performance  and  to  predictive  uncertainty  regions  for  the  location  of  an  explosive  event. 
This  project  seeks  to  determine  the  detection  performance  of  small  infrasound  arrays  and  to 
estimate  the  variances  associated  with  the  estimated  wavemmber  parameters  that  are  directly 
related  to  velocit  y  and  aziiruth.  These  characteristics  are  then  used  as  input  to  a  nonlinear 
regression  program  for  location  and  will  determine  the  overall  uncertainty  ellipse  for  locating  a 
given  event.  In  the  early  phase  of  this  project,  we  have  dev  eloped  a  methodology  for  itegrating 
w  avenuboer  estimators  from  detecting  local  arrays  into  an  overall  estimator  for  location  and  its 
associated  posterior  probability  region.  Global  capability  is  expressed  as  a  contour  map  showing 
areas  of  90%  posterior  probability  ellipses  as  a  function  of  the  expected  configuration  of  detecting 
stations  (see,  for  example,  Shumw  ay,  2000).The  expected  configuration  is  based  on  the  signal 
detection  capabilities  of  the  given  collection  of  sub-arrays  and  the  variance  properties  of  the 
maximum  likelihood  estimators  of  velocity  and  azimuth  for  each  sub-array.  The  contour  map  is 
based  on  simulations  involving  500  hypothetical  events  originating  from  each  point  on  a  5-degree 
grid  covering  the  surface  of  the  earth.  Each  event  is  detected  at  a  random  configuration  of 
International  Monitoring  Stations  (IMS). 

In  our  current  researc  h  efforts,  ve  are  focusing  on  the  problems  involved  in  evaluating  the  local 
sub-array  input  parameters  that  determine  location  uncertainty  and  the  average  areas  mentioned 
above.  The  characteristics  of  interest  are  sub-array  detection  probabilities  and  estimators  for  the 
variance  of  the  \\a\n lumber  estimates  corresponding  to  the  velocity  and  azimuth  of  the 
propagating  signal.  In  previous  reports  and  papers  (Shumw  ay  et  al,  199,  Shmw  ayl999  and 
Shumw  ay,  2000),  whave  characterized  small  array  performance  for  the  optimum  single  signal 
detector  in  terms  of  bandwidth,  signal-to-noise  ratio,  array  configuration,  and  signal  decorrelation. 
We  are  currently  applying  this  technique  to  a  number  of  even  ts  con  tained  in  the  infrasound 
database  located  at  the  Center  for  Monitoring  Research,  referred  to  in  the  sequel  as  the  CMR 
R&D  Test  Bed  Infrasound  Waveform  Library  The  database  contains  earthquakes,  explosions, 
bolides  and  missile  launches.  There  is  a  meteorite  recorded  at  six  IMS  type  stations  that  could 
provide  useful  input  parameters  for  evaluating  global  location  capabilities. 

A  second  objective  of  this  next  phase  is  investigate  the  effects  of  interfering  signals  and  (or) 
correlated  noise  on  the  conventional  F  detector,  wliih  is  based  on  the  assumptions  that  signals  are 
perfectly  correlated  and  that  noise  is  spatially  white.  We  also  consider  the  performance  of 
alternative  estimators  that  are  advertised  to  be  effective  in  a  multiple  signal  context.  Among  these 
are  the  single  and  multiple  signal  F  detectors  of  Shumw  <y  (1983),  the  estimator  of  Capon  (1969) 
and  the  MUSIC  estimator  proposed  by  Schmidt  (1979).  We  evaluate  these  alternative  estimators 
here  on  a  mixture  of  tw  o  signals  constructed  ly  adding  the  Small  Fry  nuclear  explosion,  contrived 
to  arriv  e  at  225  degrees  to  an  iherfering  signal  arriving  at  135  degrees.  Figure  1  shows  the 
primary  even  t  in  the  left  hand  column  and  the  irterfering  invent  in  the  center,  with  the  sum 
(SF+.5GS)  shown  in  the  last  column.  The  event  is  assumed  to  be  observed  at  the  three  component 
triangular  array  with  approximately  1  km  sides,  as  described  by  Shumw  av  et  al  (1998) .It  can  be 
seen  that  the  mixture  obscures  the  primary  signal  quite  effectively. 

RESEARCH  ACCOMPLISHED 


In  order  to  investigate  possible  approadres  to  evaluating  the  effect  of  contaminating  signals  such  as 
the  one  in  Figure  1,  we  consider  frequency  domain  version  of  single  and  multiple  station  signal 
detectors  with  velocity  and  aziiruth  parameterized  by  a  probe  vector  at  wavemmber  §=  (61,82)', 
i.e., 

^0)=  (1) 

where  the  array  coordinates  in  km,  relative  to  a  center  sensor  at  the  origin,  are  denoted  by 
Tj  =  (rjl,rj‘2,j  =  1, . . . ,  N-,  in  Figure  1,  there  are  N  =  3  sensors. 
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Figure  1.  Infrasound  explosive  signals  from  225°  and  135°  and  a  mixture  sampled  at 
5  Hz  at  a  triangular  array  with  1  km  sides  . 


F  or  the  classical  maximum  likelihood  solution  in  either  the  fixed  or  random  signal  case,  the 
estimated  wavemmber  vector  is  the  one  maximizing  the  beam  power,  say 

B(S)  =  ox*(0)S  x(0),  (2) 


where 


s  =  j2^n 


is  the  spectral  matrix  evaluated  at  K  frequencies  in  the  neighborhood  of  some  assumed  signal 
frequency  /o,  and  ¥*,  denotes  the  N  x  1  vector  of  discrete  transforms  of  the  observed  waveform.lt 
should  be  noted  that  (2)  is  realized  simply  online  by  filtering  the  signal  in  the  neighborhood  of  the 
center  frequency  and  delaying  and  summing  at  the  velocity  and  azimuth  matching  the  wavenunber 
vector.  A  pow  er  detector  is  then  applied  to  the  resulting  beam  to  get  an  approximation  to  (2).  All 
detectors  here  are  expressed  in  terms  of  the  spectral  matrix  S  to  provide  a  rough  comparison  to 
other  detectors  introduced  below.  Shurnw  ay  et  al  (1999)lme  also  demed  the  large-sample 
covariance  matrix  of  the  maximum  likelihood  estimator  §.  which  can  then  be  directly  incorporated 
into  the  nonlinear  least  squares  program  for  location  (See  Shurnw  ay,  2000). 

The  maximum  likelihood  detector,  derived  from  a  likelihood  ratio  test  of  no-signal  in  either  the 
random  or  fixed  signal  case  is  the  F-Statistic 


F1m  = 


ox*  (§)S  x(0)/N 
tr (S)-oaf(0)S  x(0)/N 


(N  —  1), 


which  has  an  F  distribution  with  2 K  and  2K(N  —  1)  numerator  and  denominator  degrees  of 
freedom  respectively  when  the  wavemmber  is  correct  (tr  denotes  trace).  F or  online  applications, 
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one  can  filter  the  channels  into  the  frequency  band  of  interest  and  then  use  the  ordinary  time 
domain  F-Statistic  with  2 BT  and  ‘2BT(N  —  1)  degrees  of  freedom.  The  numerator  is  proportional 
to  the  output  of  a  pow  er  detector  and  the  denominator  can  be  computed  as  the  difference  betv  een 
the  stac  ked  pw  er  and  the  beam  pav  er  in  the  numerator,  when  the  data  have  been  filtered  into  a 
frequency  band  of  width  B  Hz  and  the  sample  length  is  T  seconds.  Among  the  advantages  of  this 
detector  are  (a)  optimality  under  perfect  signal  correlation  and  spatially  white  noise,  )(b)  known 
statistical  distribution  that  does  not  depend  on  signal  and  noise  parameters  and  (c)  existence  of  an 
easily  applied  online  time  version  for  monitoring  applications.  P  oten  tial  disadantages  of  the 
simple  F  detector  are  lack  of  robustness  to  signal  correlation,  investigated  theoretically  in 
Shumw  ay  et  al  (1999)  and  interfering  signals.  These  latter  tw  o  aspects  are  assessed  here.  We  also 
study  three  alternative  estimators  that  may  have  some  potential  for  improved  estimation  in  the 
mixed  signal  case. 

The  first  is  due  to  Capon  (1969)  and  involves  maximizing  the  quadratic  form 


C(0)  = 


(5) 


over  possible  raven  miner  vectors  #.  Under  the  assumption  that  the  theoretical  spectral  matrix  of 
the  v  ectorff),  is  E,  coupled  with  the  Gaussian  assumption,  this  statistic  has  a  distribution 
proportional  to  chi-squared,  where  the  proportionality  constant  is  af(S )E_1ar(#)  (see  Capon  and 
Goodman,  1970,  for  the  exact  result).  Hence,  the  distribution  is  known  but  depends  on  nuisance 
parameters  in  the  spectral  matrix  E.  F  urthermore,  asymptotic  results  for  the  estimator# 
maximizing  (5)  are  not  available  so  there  the  inputs  required  for  assessing  global  location 
uncertainty  are  not  available,  the  statistic  is  sometimes  referred  to  as  a  maximum  likelihood 
estimator  but  the  connection  to  maximum  likelihood  is  tenuous  at  best.  If  S  were  the  maximum 
likelihood  estimator  of  the  noise  covariance  matrix  E  (it  is  not)  under  the  fixed  signal  model,  C(S) 
w  ould  be  the  maxinum  likelihood  estimator  of  the  variance  of  the  best  linear  unbiased  estimator. 

An  estimator  based  on  the  eigen  vectors  of  S  is  the  centerpiece  of  the  Multiple  Signal 
Characteristic  (MUSIC)  statistics  suggested  by  Schmidt  (1979).  A  good  summary  of  the  statistical 
properties  of  this  estimator  is  Stoica  and  Nehorai  (1989).  If  we  let  ®i, . . ..,  wjv  be  the  eigen  vectors 
of  S,  then  the  MUSIC  estimator  is  the  value  of  #  maximizing 

(6) 

where  M  <  N  is  the  assumed  number  of  signals  and  wm+  i,  ■■■,««  are  the  N  —  M  smallest  eigen 
values.  The  estimator  follows  from  the  fact  that  for  the  spatially  white  noise  model  with  correlated 
signals,  the  waven  miner  vectors  are  orthogonal  to  the  matrix  containing  the  last  N  —  M  eigen 
vectors  as  columns.  Pre-multiplying  by  af(0)  and  taking  the  inverse  leads  to  (6).  Stoica  and 
Nehorai  (1989)  have  deriv  ed  the  large  sample  coariance  matrix  for  a  linear  array. 

A  third  possibility  is  a  multiple  fixed  signal  approach  via  likelihood  ratio  tests  as  proposed  by 
Shumw  ay  (1983). This  proceeds  by  minimizing  the  squared  error 


M(S)  = 


E  ’ 

A  — 


Wj  IWj 


K 

SSE2{&)  =  \\wk  -4h)S%k  -a(ft)S2J2 

k=i 


over#  =  (#).#!>)',  Sik,  S2k,  k  =  1 .....  /\'  corresponding  to  tw  o  signals  at  wive-numbers  &i  and  #o. 
Denote  the  maximized  value  of  the  above  by  SSE 2(#)  and  the  maximized  value  with  no  second 
signal  as  SSEl(§i).  Then,  the  multiple  signal  F-statistics  for  test  the  no-second-signal  hypothesis 
is 


_  N  -  2  [SSE l(#i)  -  SSE2(§)} 
2  SSE2(§2 ) 


(7) 
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which  converges  to  an  F  distribution  with  ‘2K  and  2I\(N  —  2)  degrees  of  freedom. 

We  have  tested  the  methods  on  the  pure  signal  in  the  left-hand  column  of  Figure  1  and  on  the 
mixture  in  the  right-hand  column.  The  result  are  shown  in  Figure  2,  3  and  4  and  summarized  in 
T  able  1.  All  plots  are  constructed  using  a  center  frequency  of  .05  Hz,  corresponding  to  the  spectral 
peak  of  the  observed  data  and  K  =  13  in  the  spectral  matrix,  which  implies  a  bandwidth  of 
(13/1024)5  =  .06  Hz  so  that  the  frequency  spans  the  interval  .03  —  .09  Hz. 


Beam  Power  Capon  Noise  Power  Contours 


Figure  2.  Detectors  applied  to  pure  infrasound  explosive  signal  from  225°.  Standard 
errors  of  the  maximum  likelihood  estimators  are  in  parentheses  on  the  F- 
Statistic  contour. 


Figure  2  shows  the  wa\enuirber  plots  of  the  beam  detector  (2),  F-Statistic  (4),  Capon  detector  (5) 
and  the  MUSIC  detector  (6)  for  the  pure  signal  case.  All  statistics  give  azimuths  within  ±2°  of  the 
kno  wn  azimuth.  The  standard  deviation  of  the  maximum  likelihood  azimuth  estimator  was  8°.  F 
value  of  11.31  exceeds  the  .001  significance  point  of  the  F-Statistic  with  2(13)  =  26  and 
2(13) (3  —  1)  =  52  degrees  of  freedom.  The  peaks  in  the  beam  pow  er  and  Capon  statistics  are 
approximately  the  same  width;  the  Capon  contours  are  slightly  distorted  due  to  the  non-diagonal 
spectral  matrix  S.  The  F-Statistic  and  MUSIC  plots  promise  more  resolution,  due  to  narrow  er 
peaks. 

We  see  in  Figure  3  that  narrow  peaks  do  not  necessarily  translate  into  an  improved  ability  to 
separate  mixed  signals.  The  detectors  applied  to  the  mixture  still  focus  on  the  stronger  signal. 
Estimators  for  the  azimuth  are  biased,  ranging  from  215  to  217°,  indicating  that  the  contaminating 
signal  has  pulled  these  values  toward  its  azimuth  of  135°.  Hence,  there  is  no  appreciable  difference 
betw  een  the  bias  terms  for  the  three  methods.  The  standard  deviations  of  the  azimuth  estimators 
are  approximately  doubled  to  17°  and  the  F  value  is  reduced  to  3.55,  which  still  exceeds  the  .001 
significance  value.  It  is  clear  that  the  presence  of  the  contaminating  signal  in  this  case  increases 
the  estimated  variance,  implying  that  this  particular  array  will  have  less  influence  on  the  overall 
location. 

w  avenuhner  plots  for  the  two  multiple  signal  possibilities,  MUSIC  in  (6)  and  the  multiple  signal  F 
detector  (7)  are  shown  in  Figure  4.  The  MUSIC  estimator  (6)  in  the  bottom  panel  uses  values  of 
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Beam  Power  Capon  Noise  Power  Contours 


Figure  3.  Single  signal  detectors  applied  to  mixture  of  tw  o  pure  infrasound  explosive 
signals. 


F-Statistics:  Signal  1 


F-Statistics:  Signal  2 


Figure  4. 


Music  Eigenspectrum:  Signal  1 


Music  Eigenspectrum:  Signal  2 


Multiple  signal  detectors  applied  to  mixture  of  tw  o  infrasound  explosive 
signals. 


M  =  1  for  the  single  signal  model  and  M  =  2  for  the  tw  o-signal  lersion.  Curiously,  the  two  signal 
model  focuses  only  on  the  imbedded  signal  velocit  y  and  aziiruth,  whereas  the  single  signal  model 
peaks  for  the  stronger  signal.  The  F-Statistics  are  taken  as  F|  (&\ )  in  the  left  hand  column  and 
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To  (Hi.  #2)  in  the  second  column,  where  we  assume  the  biased  value  for  B\  in  the  computation  of  the 
statistic  (7).  The  F-Statistic  for  the  second  signal,  given  the  first  signal  at  217°  is  2.63,  which  still 
exceeds  the  .01  false  alarm  level  of  2.55  with  2(13)  =  26  and  2(13) (3  —  1)  =  26  degrees  of  freedom 
in  the  numerator  and  denominator. 

A  convenient  summary  of  the  results  is  given  in  Table  1  and  we  note  the  relative  robustness  of  the 
single  station  F  detector.  Not  also  that  despite  a  bias  of  about  8°  in  the  maximum  likelihood 
estimator,  the  increase  in  standard  deviation  still  puts  the  estimated  value  within  half  a  standard 
deviation  of  the  ground  truth  azimuth  of  225°. 


T  ablel:  Estimation  by  Different  Methods  for  Single  and  Mixed  Events  Using  Single  and  Multiple 

Signal  Detectors 


Event 

SF225  Only 

Mix 

Detectors 

SF225 

SF225 

GS135 

Single  F 

225(8) 

217(17) 

- 

F- Values 

11.31 

3.55 

- 

Capon 

227 

215 

- 

Multiple  F 

225 

217 

134 

F- Values 

11.31 

3.55 

2.63 

MUSIC 

223 

216 

131 

CONCLUSIONS  AND  RECOMMENDA  TIONS 

We  conclude  that  the  performance  of  the  likelihood  based  statistics  under  the  assumptions  of 
spatially  white  noise  and  possible  signal  decorrelation  are  relative  robust  to  departures  from  these 
assumptions.  We  have  also  examined  the  Capon  and  MUSIC  detectors  as  possible  competitors  to 
the  F-Statistic  and  have  note  that  they  do  not  improve  performance  in  the  mixed  signal  case  and 
do  have  convenient  expressions  for  the  variances  of  the  estimated  wavemmber  parameters.  We 
note  that  the  likelihood  based  approaches  produce  large  sample  variances  for  the  estimated 
w  avenuber  parameters  (velocity  and  azimith)  that  are  needed  for  input  into  the  nonlinear 
location  simulations. 

A  reasonable  number  of  calibration  even  ts  are  needed  to  giv:  realistic  inputs  for  the  previously 
developed  softw  are  that  produce  area  uncertairty  contours  for  the  location  of  events  detected  by 
random  configurations  of  detecting  stations  in  the  IMS  netw  ork.for  these  inputs,  we  are  currently 
downloading  a  number  of  events  con  tained  in  the  CMR  R&D  Tst  Bed  Infrasound  Waveform 
Library.  Events  available  are  a  California  earthquak,  several  gas  pipe  explosions,  a  Titan  IV  B 
launch  and  some  bolides. 
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